Effect of heat flux on difl'erential rotation in turbulent convection 
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We studied the effect of the turbulent heat flux on the Reynolds stresses in a rotating turbulent 
convection. To this end we solved a coupled system of dynamical equations which includes the 
equations for the Reynolds stresses, the entropy fluctuations and the turbulent heat flux. We used 
a spectral r approximation in order to close the system of dynamical equations. We found that 
the ratio of the contributions to the Reynolds stresses caused by the turbulent heat flux and the 
anisotropic eddy viscosity is of the order of ~ 10 {Lp/lo)^ , where lo is the maximum scale of turbulent 
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differential rotation and should be taken into account in the theories of the differential rotation of 
the Sun, stars and planets. In particular, we demonstrated that this effect may cause the differential 
rotation which is comparable with the typical solar differential rotation. 



PACS numbers: 47.27.-1; 47.32.-y 



I. INTRODUCTION 



Solar and stellar magnetic fields are believed to origi- 
nate in a dynamo, driven by the joint action of the mean 
hydrodynamic helicity of turbulent convection and dif- 
ferential rotation (see, e.g., [E0,|1,0,H|, and references 
therein) . It was suggested in [g] that the differential rota- 
tion of the Sun is caused by an anisotropic eddy viscosity 
which was described phenomenologically in 0, H, IM- 
Theory of the differential rotation based on the idea of 
the anisotropic eddy viscosity was developed in a number 
of papers (see, e.g., 0,0,13, and references therein). 
However, there is an additional effect which can strongly 
modify the differential rotation. In particular, the direct 
effect of the turbulent heat flux on the Reynolds stresses 
in a rotating turbulent convection is crucial for formation 
of the differential rotation. 

The effect of rotation on a hydrodynamic turbulence 
was studied in numerous papers (see, e.g., 0,0|). How- 
ever, a relation to the turbulent convection was made in 
previous theories of the differential rotation only phe- 
nomenologically, using the equation 



cx gTo{u'^s') 



(1) 



which follows from the mixing-length theory. Here {u'^s') 
is the vertical turbulent heat flux, u' and s' are fluctua- 
tions of fluid velocity and entropy, g is the acceleration 
of gravity and tq is the characteristic correlation time 
of turbulent velocity field. Equation implies that 
the vertical turbulent heat flux plays a role of a stirring 
force for the turbulence. However, a more sophisticated 
approach implies a solution of a coupled system of dy- 
namical equations which includes the equations for the 
Reynolds stresses (u^u'), the turbulent heat flux (s'u^) 
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and the entropy fluctuations (s'^) in a rotating turbulent 
convection. The latter has not been taken into account 
in the previous theories of the differential rotation. 

The goal of this study is to analyze the effect of the 
turbulent heat flux on the Reynolds stresses in a rotating 
turbulent convection and on formation of the differential 
rotation. We demonstrated that this effect is crucial for 
the formation of the differential rotation, and it should be 
taken into account in theories of the differential rotation 
of the Sun, stars and planets. In particular, we found that 
the ratio of the contributions to the Reynolds stresses 
caused by the turbulent heat flux and the anisotropic 
eddy viscosity is of the order of '--^ 10(Lp/Zo)^ for tq < 1, 
where is the maximum scale of turbulent motions (the 
integral scale of turbulence), is the rotation rate, Lp is 
the fluid density variation scale, i.e., {'Vpo)/po = — L^^e, 
Pq is the fluid density and e is the unit vector in the di- 
rection of the fluid density inhomogeneity. The turbulent 
heat flux contribution to the Reynolds stresses changes 
its sign when the direction of the vertical turbulent heat 
flux changes. This is the key difference from previous 
theories of the differential rotation. The effect of the tur- 
bulent heat flux on the Reynolds stresses in a turbulent 
convection may cause the differential rotation which is 
comparable with the typical differential rotation of the 
Sun. The data of the solar differential rotation is ob- 
tained from surface observations of the solar angular ve- 
locity (see, e.g., Bind from helioseismology based 
on measurements of the frequency of p-mode oscillations 
(see, e.g., E [13 El 111). 

The mechanism of the differential rotation that is as- 
sociated with the effect of the turbulent heat flux on 
Reynolds stresses in a rotating turbulent convection is 
as follows. Let us split the total rotation of fluid into 
a constant component f2 (uniform rotation) and the dif- 
ferential rotation SQ. The uniform rotation causes the 
counter-rotation turbulent heat flux (i.e., the toroidal 
turbulent heat flux, (s' u'^), directed opposite to the back- 
ground rotation fi). Therefore, there is a correlation of 
fluctuations of the entropy s' and the toroidal compo- 
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nent of tfie velocity, u'^. Here r,9,ip are the spherical 
coordinates. 

The counter-rotation turbulent heat flux is similar 
to the counter-wind turbulent heat flux (in the direc- 
tion opposite to the mean wind) which is well-known 
in the atmospheric physics. The counter-rotation tur- 
bulent heat flux arises by the following reason. In tur- 
bulent convection an ascending fluid element has larger 
temperature then that of surrounding fluid and smaller 
toroidal fluid velocity, while a descending fluid element 
has smaller temperature and larger toroidal fluid veloc- 
ity. This causes the turbulent heat flux in the direction 
opposite to the toroidal mean fluid flow (i.e., opposite 
to rotation). The counter-rotation turbulent heat flux is 
determined by Eq. in Section III. 

The entropy fluctuations cause fluctuations of the 
buoyancy force, and this results in increase fluctuations 
of the vertical and meridional components of the ve- 
locity which are correlated with the fluctuations of the 
toroidal component of the velocity. These produce the 
off-diagonal components of the Reynolds stress tensor, 
{u'j.u'^) and {u'gu'^), and create the toroidal component 
of the effective force which causes formation of the dif- 
ferential rotation 6fl in turbulent convection. 

This paper is organized as follows. In Section II we for- 
mulated the governing equations, the assumptions, the 
procedure of the derivation, and described the effect of 
the turbulent heat flux on the Reynolds stresses. In Sec- 
tion III we developed the theory of the differential rota- 
tion based on this effect. In Section IV we made esti- 
mates for the solar differential rotation. In Appendixes 
A and B we performed a detailed derivation of the effect 
of the turbulent heat flux on the Reynolds stresses in the 
rotating turbulent convection. 

II. EFFECT OF THE TURBULENT HEAT 
FLUX ON THE REYNOLDS STRESSES 

In order to study the effect of the turbulent heat flux 
on the Reynolds stresses we considered turbulent con- 
vection with large Rayleigh and Reynolds numbers. We 
employed a mean field approach whereby the velocity, 
pressure and entropy are separated into the mean and 
fluctuating parts, where the fluctuating parts have zero 
mean values. The large-scale fluid motions are deter- 
mined by the mean-field equations, which follow from 
the momentum and entropy equations for instantaneous 
fields by averaging over an ensemble of fluctuations. The 
mean-field equations are given by 

(^^ + U- [/, = -V, -g5 + 2Ux 

- {V,+L;'e,){u',u'^)+U^) , (2) 

- -iv-F„(U,5) , (3) 



where Eq. ^ is written in the reference frame uniformly 
rotating with the angular velocity CI. Here the mean 
fields U, P, T and S are the fluid velocity, pressure, 
temperature and entropy, respectively, po (U) is the 
mean molecular viscous force, Fk(U, S) is the mean heat 
flux that is associated with the molecular thermal con- 
ductivity. The mean fluid velocity U for a low Mach 
number flows satisfles the equation div (po U) = 0. Equa- 
tions (0) and (PJ are written in the anelastic approxima- 
tion, which is a combination of the Boussinesq approxi- 
mation and the condition div (po U) = 0. The variables 
with the subscript "0" correspond to the hydrostatic 
nearly isentropic basic reference state, i.e., VPq = PoS 
and g • [(7Po)^^ VPq — Po ^ Vpo] ~ 0, where 7 is the ratio 
of specific heats. The turbulent convection is regarded as 
a small deviation from a well-mixed adiabatic reference 
state. 

In order to get a closed system of the mean-field 
equations we have to determine the dependencies of 
the Reynolds stresses (u-(i, x) u'j{t, x)) and the turbulent 
heat fiux (s'(t, x) u[{t, x)) on the mean fields. To this end 
we used equations for fluctuations of velocity and entropy 
in a rotating turbulent convection, which are obtained 
by subtracting Eqs. 10) and (PJ for the mean fields from 
the corresponding equations for the instantaneous fields. 
The equations for fiuctuations of velocity and entropy are 
given by 

-g.s' + 2u' X n + U^ , (4) 
% = -^iu'-e)-{U-V)s' + S\ (5) 

where = ((u' • V)u') - (u' • V)u' + f^(u') and 
5^ = ((u • V)s) - (u • V)s - (l/To) V • F,(u',s') are 
the nonlinear terms which include the molecular dissipa- 
tive terms, fi^ = — g • VS* is the Brunt- Vaisala frequency, 
p' are fiuctuations of fiuid pressure and the fiuid velocity 
fiuctuations u' satisfy the equation div (po u') = 0. 

To study the rotating turbulent convection we per- 
formed the derivations which include the following steps: 
(i) using new variables for fluctuations of velocity v — 
^/pou' and entropy s — y^s'; (ii) derivation of the 
equations for the second moments M^-^^^ of the velocity 
fluctuations {vi Vj), the entropy fluctuations (s^) and the 
turbulent heat flux {vi s) in the k space; (iii) application 
of the spectral closure (see Eq. below) and solution 
of the derived second-moment equations in the k space; 
(iv) returning to the physical space to obtain formulae 
for the Reynolds stresses and the turbulent heat fluxes 
as the functions of the rotation rate H, (see for details. 
Appendix A). 

The second-moment equations include the first-order 
spatial differential operators A/" applied to the third-order 
moments M^^^^\ A problem arises how to close the sys- 
tem, i.e., how to express the set of the third-order terms 
A/'M^^^^^ through the lower moments M*^^^) (see, e.g.. 
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[20l I21L 12^'). Various approximate methods have been 
proposed in order to solve it. A widely used spectral r 
approximation ( US El IH IH IH 113 ) postulates that 
the deviations of the third-moment terms, AfM^-^-^-^^li), 
from the contributions to these terms afforded by the 
background turbulent convection, AAMq^^^'' (k) , are ex- 
pressed through the similar deviations of the second mo- 
ments, M(")(k) - M^"\k): 



Af(^^)(k) - Af(^'''(k) 



(6) 



where T(fc) is the characteristic relaxation time, which 
can be identified with the correlation time of the turbu- 
lent velocity field. The background turbulent convection 
(which corresponds to a nonrotating and shear-free tur- 
bulent fluid flow) , is determined by the budget equations 
and the general structure of the moments is obtained 
by symmetry reasoning. The above procedure (see Ap- 
pendix A) yields formulae for the Reynolds stresses and 
the turbulent heat flux in the rotating turbulent convec- 
tion. In particular, this allowed us to determine the effect 
of the turbulent heat flux on the Reynolds stresses. 

The differential rotation in the axisymmetric fluid flow 
is determined by linearized Eq. (01 for the toroidal com- 
ponent U^{r, 9) = r sinOSfl of the mean velocity: 



Po ■ 



du^ 15,3, I 

ot r'^ or r sm' 

+2pQ (Uxn)^ , 



d 
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(7) 



where the tensor tr, 



(viVj) is determined by the 



Reynolds stress tensor. In particular, 



eT ^"i '"3) = Po 



T dr\r 



. F I u 
rep ~^ ru:> ' 



_ V e /„, „, \ ,, sm9 d f 



ae^ = -ej {v, vj) = po — 



(8) 



(9) 



where e'', and e*^ are the unit vectors along the ra- 
dial, meridional and toroidal directions of the spherical 
coordinates r, 0, ip. There are three contributions to the 
tensor Uij — —{viVj). In particular, the first term in the 
right hand side of Eqs. (|HJl and © describes the isotropic 
turbulent viscosity oc i/^, the second term in Eqs. 
and lO determines the contribution cr^ to the Reynolds 
stresses caused by the turbulent heat flux, and the third 
term in Eqs. ||SJ) and © determines the contribution cr" 
to the Reynolds stresses caused by the anisotropy of tur- 
bulence due to the nonuniform fluid density and uniform 
rotation. 

In Eq. |(7J) we neglected small molecular viscosity term 
and we took into account that in the axisymmetric fluid 
flow 9(U, S) /dip = and Si^/il is a small parameter. We 



assumed that the toroidal component of the mean veloc- 
ity is much larger than the poloidal component. This 
is typical for the solar and stellar convective zones. We 
also took into account that the fluid density is nonuni- 
form in the radial direction. The first two terms in the 
right hand side of Eq. Q are the (p component of the 
divergence of the tensor (Jij written for the axisymmetric 
fluid flow in spherical coordinates. This is a standard 
form of the ip component of the divergence of a tensor in 
spherical coordinates (see, e.g., psf). 

Let us discuss the contribution, <t^, to the Reynolds 
stress tensor caused by the turbulent heat flux. The 
components of this tensor, cr^ = —e^ el {viVj)-^ and 
— — eJ {vi Vj)^ , were determined in Appendix A. 
They are given by 



^poT^gF^n sin6l[$i(u;) + cos^ 6* $2(w)] , 



- PqTq g F^Q sin^ 9 cos 6* $2 (w) , 



(10) 



(11) 



where the tensor {vi vj)^ determines the contribution to 
the Reynolds stresses which vanishes when {s' u'^) 
(see Eq. (|A24|l in Appendix A). Here = {s' u'^ is the 
vertical turbulent heat flux, the parameter lo = 8f2ro, 
To = Zo/ ""0 is the characteristic correlation time of turbu- 
lent velocity field, and uq is the characteristic turbulent 
velocity in the maximum scale of turbulent motions ^o- 
The functions and ^2{'^) are shown in Fig. 1. The 

formulae for the functions and <I'2(^) are given by 

Eqs. ljA28|l and (|A29|I in Appendix A. The asymptotic 
formulas for cr^ and for a slow rotation, 8 tq il <C 1, 
are given by 

Po Tq gF^n sin 9 , 
4 Po T-Q sin^ 9 cos 9 , 

and for 8 tq 3> 1 they are given by 
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Po g sin 9 , 
Po g F^ sin^ 9 cos 9 



The contribution to the Reynolds stresses, cr", in 
Eqs. (|H1) and lO is caused by the anisotropy of turbu- 
lence due to the inhomogeneous fluid density and uni- 
form rotation. The tensor cr" determines the anisotropic 
eddy viscosity tensor with the nonzero off-diagonal com- 
ponents which are given by 



Tip 



15L 



^{pouD (nro) sm9[ai{n)~a2in) cos"" 9] 



15L2 



(poug) (r^To) 02(^2) sin^ 61 



COS( 



(12) 
(13) 



(see, e.g., ^3); where ai{n) ^ —1 and 02(^1) ~ 2(17 to)^ 
for nro < 1, and ai(ri) OKHtq)"^] and a2{fl) ~ 



4 




FIG. 1: The functions $i(r2ro) (solid) and -l>2(nTo) 
(dashed). 



— (17ro)~^ for Htq ^ 1. The contribution to the 
Reynolds stresses, cr", due to the anisotropic eddy vis- 
cosity is smaller than that of cr^ due to the turbulent 
heat flux. In particular, the ratio \a^/a'^\ ~ 10 (Lp/lo)^ 
for f^To < 1. Note that for fast rotation rates {^Itq ^ 1) 
the validity of Eqs. ifT^ and (fT!^ is questionable, because 
the quasi-linear approximation used in [lo| does not valid 
for firo > 1. 

For large Rayleigh numbers the contributions of 
anisotropics decrease in small scales. On the other hand, 
the rotation introduces an anisotropy in the turbulent 
convection. This causes non-zero off-diagonal compo- 
nents of the Reynolds stress tensor (see Eqs. (fTn|) - ([T^ ) 
and results in the redistribution of the turbulent heat flux 
on the surface of the rotating body (see Eq. l)A30|l ). Note 
that the main contribution to the tensor aij is at the max- 
imum scale of turbulent motions. Therefore, the contri- 
butions to the tensor aij which depend on the Reynolds 
number are negligibly small. In the present study we 
investigated the large-scale effects (the differential rota- 
tion), and the influence of the molecular viscosity and 
molecular thermal diffusivity on the large-scale dynam- 
ics is very small in comparison with that of the eddy 
viscosity and turbulent thermal diffusivity. Therefore, 
the contributions to the tensor ct^ which depend on the 
molecular Prandtl number are also negligibly small. 



III. DIFFERENTIAL ROTATION 

In the present study for simplicity we have taken into 
account only the effect of cr^ on the differential rota- 
tion. Let us neglect the toroidal component of the Cori- 
olis force in Eq. (|7|). This is valid when the poloidal 
components of the mean fluid velocity Ug and Ur are 
much smaller than uolo/Lp. This condition implies that 
the last term in the right hand side of Eq. (Q is much 
smaller than other terms. Here we took into account that 



cr^ cx pq and the fluid density stratification length Lp is 
much smaller than the solar radius Rq. Therefore, we 
neglected the effect of the meridional circulations on the 
differential rotation, which was studied in [l^ . among 
the others. For simplicity we also did not take into ac- 
count the dependence of the turbulent viscosity on the 
rate of rotation. After these simplifications, Eq. {Tj) in 
dimensionless form reads 



where 



d_ 

dt 



Mix) 

W(r) 
X(r) 



M{X) 



W{r) 



sn 



(14) 



AX 



dx^ ' "'ax ' 

Po{r) dr \ dr 

/o[(3-rA)$i(co)-2$2(t^) 
+ (13-rA)$2(t^)X2] , 



^0 = g F^:/Qv^, A — R^fLp and X = cos9. Here 
length is measured in units of the solar radius Rq and 
time is measured in units of Rq/v^ based on the solar 
radius and the turbulent viscosity i^^ . 

Solution of Eq. p4|l we seek for in the form: 



sn . 
77 = ^ 



1 



C: 



3/2 



n=0 



(X) VmMr)QmMt) 
m— 3 



(15) 



where the function Cn^'^{X) satisfies the equation for the 
ultra-spherical polynomials: 



[M{X) - n{n + 3)] C^/^{X) ^ . 



(16) 



The functions Vm.niT) in Eq. H15|l are determined by 
equation of the eigenvalue problem 



[Cnir) - 7m]K^,„(r•) = 0, (17) 



with 



^ , , £ 4 d 2A A2 n(n + 3) 

^n{r) = T2+-:r + -A 

ar'^ r dr r 4 r- 

The constant A in Eq. H15|l is determined from the con- 
servation law for the total angular momentum Lq = 
/{rx[(ri -I- 5Vl)xv\} pq{v) dv of the rotating body (e.g., 
the Sun). The functions Qm,2n{t) in Eq. p4|l are deter- 
mined in Appendix B. They are given by Eqs. (|B4p - (|B6p . 

Equation (|17|l coincides with the equation for the Ke- 
pler problem in quantum mechanics (the hydrogen atom 
in a spherically symmetric potential, see, e.g., |29j|). The 
solution of Eq. H17|l is given by 

V„,,n{r) = r" exp(-r A/m)i?(a; b; 2rA/m) , (18) 

where F{a] b; y) is the confluent hypergeometric function 
with a = n — m -\- 2 and b = 2(n -f 2). Here we assumed 
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for simplicity that Lp is independent of the radius. The 
characteristic spatial scale of the mean fields variations 
is of the order of the solar radius i?© in the main part 
of the convective zone (except for its boundaries). On 
the other hand, the fluid density in the solar convective 
zone changes very strongly (in 6 - 7 orders of magnitude). 
Therefore, in the main part of the solar convective zone 
rA 1, and Eq. (fTS|l for the eigenfunctions Vm,n{f') 
reduces to 



where the eigenvalues 7„ 

^711 



' ^ exp(— r A/to) , 
are given by 
A2(m2 -4) 



Am? 



(19) 



(20) 



with the integer numbers to > 3. Here 7™ is measured 
in units of v^/R^. The constant A in Eq. H19|) is deter- 
mined from the condition „ (r) dr = 1 , where 
fb = Rb/ Rq , Rb is the radius of the bottom of the con- 
vective zone, and r = R/Rq is the dimensionless radius 
measured in units of the solar radius Rq. Therefore, the 
differential rotation caused by the effect of the turbulent 
heat flux on the Reynolds stresses in a turbulent convec- 
tion is determined by the following equation: 



sn 



Lp 
Ro 



Po{Rb 



/ R Y"-' 


'po{Rb)' 


\RcJ 


[ PoiR) \ 



Pq{Rq) 

-_2 

[/l,m(^) 



+/2,m(w) COS^O] 



(21) 



where 

/l,m(t^) 



K{m) 



l>l(co)-3if,(TO)^|.2(^) 



K{m) $2(^) 



10 TO 



fl + ./2/5 



1 - (Rb/Rer 



Rq \m - 

K{m) = 8TOV[3(m2 - 4)(to -|- 2)], $i(w) = $i(tj) + 

(13/10)<i>2(w) and ^2{i^) = + 5$2(^), the param- 

eter K^:(ra) is given by Eq. HB7() in Appendix B. In the 
next section we use Eq. (I21II in order to estimate the solar 
differential rotation. 

IV. DISCUSSION 



little difference to these estimates. In particular, at depth 
of the convective zone, H ~ 10^" cm measured from the 
top (i.e., at i? = O.SSi?©), the parameters are: the max- 
imum scale of turbulent motions Iq ~ 5.5 x 10^ cm; the 
characteristic turbulent velocity in the maximum scale 
of turbulent motions uo ~ 5.4 x 10"^ cm s^^; the tur- 
bulent viscosity ~ 10^'^ cm^ s~^; the fluid density 
po ~ 7.6 X 10~^ g cm~'^ and the fluid density stratifi- 
cation length Lp ~ 10^° cm. Thus, Eq. yields the 
following estimates for the solar differential rotation 



d_ 
dr 
d_ 
89 



6n 



6n 



200 (1 + 0.4 cos^ 



nHz 

Rrr, 



140 sin(26') 



nHz 

Rp) 



(22) 
(23) 



These estimates are in agreement with the data obtained 
from surface observations of the solar angular velocity 
Ullil and from hehoseismology E IH El Ei . There- 
fore, the effect of the turbulent heat flux on the Reynolds 
stresses in a turbulent convection is crucial for the for- 
mation of the differential rotation and should be taken 
into account in theories of the differential rotations of the 
Sun and solar-like stars. 

The mechanism of the differential rotation due to the 
effect of the turbulent heat flux on the Reynolds stresses, 
is related to the counter-rotation turbulent heat flux in 
turbulent convection. This flux reads: 



( arctan 
2 1 

. v ^ 



-f ln(l+tj2) 



(24) 



[see Eq. (jA30|l ]. Therefore, the entropy fluctuations cor- 
relate with the toroidal component of the velocity. The 
entropy fluctuations result in fluctuations of the buoy- 
ancy force, that increases fluctuations of the poloidal 
components of the velocity (which are correlated with 
the fluctuations of the toroidal component of the veloc- 
ity). These produce the off-diagonal components of the 
Reynolds stress tensor which cause the formation of the 
differential rotation. 
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The effect of the turbulent heat flux on the Reynolds 
stresses in a turbulent convection may cause the differ- 
ential rotation comparable with the typical differential 
rotation of the Sun. Indeed, let us use estimates of gov- 
erning parameters taken from models of the solar convec- 
tive zone, e.g., [sol Isij. More modern treatments make 



APPENDIX A: THE REYNOLDS STRESSES IN 
ROTATING TURBULENT CONVECTION 

We use a mean field approach whereby the velocity, 
pressure and entropy are separated into the mean and 
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fluctuating parts, where the equations in the new vari- 
ables for fluctuations of velocity v — JjJ^ u' and entropy 
s = yfpos' follow from Eqs. 0) and ©: 



1 av(x,<) 



dt 



9s(x, t) 
dt 



= -V 



1 



G*^ V - g s 



2vx r2-(v V)U 
(Al) 



Pa 

Fm + , 



^^(ve)-G^s + s^ , (A2) 



where G'^ = (l/2)divU + U-V, and are the 
nonlinear terms which include the molecular viscous and 
dissipative terms, p' are fluctuations of fluid pressure. 
The fluid velocity fluctuations v satisfy the equation 
V-v = (l/2ip)(ve), where (Vpo)/po = -^p^e. Equa- 
tions IjAip and (|A2I) are written in the anelastic approxi- 
mation. The variables with the subscript "0" correspond 
to the hydrostatic nearly isentropic basic reference state. 
The turbulent convection is regarded as a small deviation 
from a well-mixed adiabatic reference state. 

Let us derive equations for the second-order moments. 
For this purpose we rewrite the momentum equation and 
the entropy equation in a Fourier space. In particular. 



dt 

ds{k) 
dt 



[DIM + /L(k)]«,„(k) -t- ge,„ P„„(k) s{k) 
+«f (k) , (A3) 



= -G^(k)s(k)+s 



N 



(A4) 



To derive Eq. (jA3(l we multiplied the momentum equa- 
tion written in k-space by Pij (k) = Sij — kij in order tc 
exclude the pressure term from the equation of motion 
Here 



/^(k) = 2fc„V,C/„- V,C/,-G^(k)5,, , 
G'^(k) = idivU + i(U-k) , 

and 5ij is the Kronecker tensor, kij = kikj/k'^ and Sijk is 
the Levi-Civita tensor. Using Eqs. ljA3(l - (|A4(l we derive 
equations for the following correlation functions: 

/ii(k) = L{v^;vj) , FiCk) ^ L{s;vi) , 
e.(k) = L{s;s) , 



where 



L{a; c) 



;a(t,k + K/2)c(t,-k + K/2)) 
X exp (iK-R) dK , 



R and K correspond to the large scales, and r and k 
to the small ones. Hereafter we omitted argument t and 



R in the correlation functions. The equations for these 
correlation functions are given by 



g^-(k) 

dt 
dF,{k) 

dt 

de{k) 
dt 

where 



u 

ijmn 



)fmn+M' +J^h, (A5) 



ijmn J J mn 



{Ig, + D2^)F^ + 5e™P™(fc)e(k) +UF, , 



= -div[ue(k)] +A/'e 



(A6) 
(A7) 



ijmn 



(k2)4 



mp 



'dk„ 



-5^m5 jn [divU + U-V] , 
DZ,{k,)5,n+D%{k2)5,n 

'-■jmp ^in) 



2 ^qkpq(^€ 



qf^pq\c-imp ^jn 



(A8) 



(A9) 



= ge„[P„„(k)F,(k)+P,„(k)P,(~k)],(A10) 

and ki = k + K/2, ka = -k + K/2. Note that the 
correlation functions fij, Fi and Q are proportional to 
the fluid density /9o(R)- Here N fij, MFi and MQ are 
the terms which are related to the third-order moments 
appearing due to the nonlinear terms. In particular. 



P„„(ki)<(ki)z;,(k2)) 



+ {v,{k,)P,MvZ{^2)) , 



(kiK(k2)) + (s(ki)P,„,(k2)«^(k2)), 



NF, = 

NQ = (S^(ki)s(k2)) + (S(ki)s^(k2)) . 

coincides with that derived 



When div U = 0, Eq. 
in 32]. 

The equations for the second-order moments contain 
high-order moments and a closure problem arises (see, 
e-g-, [13> 1^ 113 )■ We apply the spectral r approxi- 
mation (or the third-order closure procedure, see, e.g., 
[23, S M m m 123). The spectral r approxima- 
tion postulates that the deviations of the third-order- 
moment terms, M fij{k), from the contributions to these 
terms afforded by the background turbulent convection, 
jV'/j'°-' (k) , are expressed through the similar deviations 
of the second moments, fijik) — flj\k), i.e., 



J^f,,{k)-Ufg\k) 



/,,-(k)~/ff)(k) 
r(fc) 



(All) 



and similarly for other tensors, where Affij = M fij + 
M[.{F^=^) and A^P^ = NF, + ge„P„(fc)e"=o, the su- 
perscript (0) corresponds to the background turbulent 
convection (i.e., a nonrotating turbulent convection with 
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VilTj = 0), r(fc) is the characteristic relaxation time of 
the statistical monients. The quantities F^=o and 0^^=*^ 
are for a nonrotating turbulent convection with nonzero 
spatial derivatives of the mean velocity. Note that we 
applied the r-approximation (|A11|) only to study the de- 
viations from the background turbulent convection which 
are caused by the spatial derivatives of the mean veloc- 
ity and a nonzero rotation. The background turbulent 
convection is assumed to be known (see below). 

Solution of Eqs. (|A5(I - (|A7(I after applying the spectral 
T approximation reads 



/y(k) 

e(k) 



where 



4(k) 
/^(k) 
M^(k) 



[i-T(u-v)]eW(k) , 



ijmn\ / mnpqUpq ' J pq\ f 

ge„,{P,„(k)[F,(k)-Fj^="(k)] 



(A12) 
(A13) 
(A14) 

(A15) 
(A16) 



background turbulent convection is assumed to be non- 
helical. 

Equations ifHa)! . l|IT6)l and (|IT7jl can be rewritten in 
the form 



P,(k) = po F., [e,„P™(k) + ^ (exk),] , 



where P* 



1-1- 



(A23) 



2i:M^f) + 2V^(Bi - 2) M^f'] , (A24) 



M-(k) = po r(k) g '-f^ [M^f - 2^mJ)] 



(A25) 



(exk)^ Cttt, -Pj>77,(k) + (exk)j _Pjm(k) , 

Pi7n (k) -fjn (k) , 

(exk)i (exk)j , 



P,„(k)[F,(-k) - ^;^'-"(-k)]}. (A17) and we used the identities: 



Here D^-'^lJl) is the inverse of Dij{ft) 



,jmni^) is the inverse of Ajmn 
and 



rDg and 



(Jl) = 5imSjn — T D. 



ijmn ' 



^„.i(rj) = x(^) (% + ^ e^,m fcm + V'' hj) , (A18) 



-^3 {^ipjn^jn 



-\-B^ Sipm^jqn^pq -^6 (y^ipm^jpn H~ "^jpn^tpm)] i (A19) 



and fc, = fc,/fc, = 1/(1 + V-'), = 2T(fc) (k • n)/k, 

Si = l + X(2^), = Bi + 2-4x(V), B3 = 2V^x(2V'), 
S4 - 2x(V') - Si, B5 = 2 - Bi and = 2^- [x(V') - 
x(2-(/')]. For derivation of Eqs. HA12I) and (|A15l) - (IA16p we 
used a procedure described in Appendix B in 33] . 

For the integration in k-space of the second moments 
/y (k), Fi(k) and 9(k) we have to specify a model for 
the background turbulent convection (i.e., a nonrotating 
turbulent convection with WiUj — 0). Here we used the 
following model of the background turbulent convection: 



i^(°)(k) 

e(°)(k) 



Po{{nr)P^A^)W{k) , (A20) 

3po(s'<)e,„P,™(k)VK(fc) , (A21) 
2po{{sr)W{k), (A22) 



where W{k) = E{k)/8'Kp, T{k) = 2rof(/c), E{k) = 
~df{k)/dk, f{k) = (fc/fco)i"«, 1 < g < 3 is the ex- 
ponent of the kinetic energy spectrum (g = 5/3 for Kol- 
mogorov spectrum), k^ — l/lo, and is the maximum 
scale of turbulent motions, tq = Iq/uq and uq is the char- 
acteristic turbulent velocity in the scale Iq. Motion in the 



(a) 



.[{B,-B^)M,^ 
-2P3(Af,ff-Mi^))], 



(b) 



r(a) 



In order to integrate over the angles in k-space we used 
the following identities: 



(a) = 



fci, sin0 ,^ , -j- ^ - 

d0dip = AiSij + Ai uj.ij , 



1 -|- a cos^ 

ki 



sm( 



(A26) 



1 + a cos^ t 

~\-5jm^in + Sjni^im + '^nm'^ij) , (A27) 

(see for details, [Slsl)- 

where ojzj = ^l^ r^j /il^, 



1 '^tjmn 



Ai{a) 
A2{a) 



2tt 
a 

a 



(« + l) 



arctan(y^) 



^^^g^ arctan(V^) _3 



(a + 1) 



2 arctan(-y/a) 5a 



a 3 
(72(a) = ^(a)-7Ai(a)-H35Ci(a) , 
C3(a) = ^(a)-5Ci(a) . 

Equation IIA24p after the integration in k-space allows 
us to determine the contributions of the turbulent heat 
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flux to the Reynolds stress tensor in turbulent convec- 
tion. In particular, the components cr^ = — ej e[ 



"-"^ el f[. are given by Eqs. JTU)) and ((TT|l . 



and af^ = 
respectively, where 



"3 J •'U 



$i(w) = 2*i(w) + *2(t^/2) , 



(A28) 
(A29) 



arctano;,^ Scj^ 
[jj 6 



+2ujS{uj) 



arctan ui , 8cj 
(1 + cj2) + _ - 5 



cij = 8x0^2 and S{lS) — J^" [arctan y/y] dy. Here we took 
into account that = + A^^^(a)) - cf ^(w) - 

(w) and ^'2(t^) = -Cf ^(cj) - 3Cf^(cj), where the 
functions and C^^^(w) are given by 



Jo 

CiP\cj) = (6Aa;P+i) C y^Ck{y^)dy, 



(see for details, |33,IM|). For derivation of Eqs. l(Tn|l and 
we used the following identities: 

pV pS fl4-(c) _ . . cp 9 

(ef ej - ef e^) ej;,a)„ Jy„m = {C2 + SCa) sin^ 6* cos 6*, 



7- 7- 7- ^ 

^pgj 



+ cos2 6i(C'2 + 3(73)] 



Jijmn — (exa;)p[Ci + C3 



Equation for the turbulent heat flux follows from 
Eq. ljA23|l after integration in k-space: 



El 

16 



: [A[^^ (x) + 4^ (x)] sin e'^ + [2AY' (a 



+4"' (x) sin' ^] e'' + ^ ^""^ (x) sm{2e) e' 



(A30) 



where e*", and e*^ are the unit vectors along the radial, 
meridional and toroidal directions, respectively. The first 
term in Eq. ljA30|) describes the counter-rotation turbu- 
lent heat flux, which is given by l|24|l . 



APPENDIX B: THE FUNCTIONS Qm,2n{t) 

Equation for the functions Qm,2nit) follows from 
Eqs. (djl-lCTI). In particular, 



dQmfijt) 

dt 



dQma{t) 
dt 



dQm,2n{t) 
Jt 



— |7m| Qm.O + ^0 



Rr. 



i^i(m)$3(w) 



-\lrn\Qm.2~ —Io^2{lo) 
15 

-13if2(m) , 

~ I'Tjti I Qm.2n 5 



(Bl) 



— Ki (to) 



(B2) 
(B3) 



with the integer numbers 71 > 1 in Eqs. ljB3|l . where 
= -$i(cj) - (3/2)$2(^), ^4(^ ) = 3$i(w) + 
(35/2)$2(c^), Ki{m) = rVPo(r)F™,2„(r)dr and 
K2{m) = J^^ \/po{r) Vm,2n{r) dr. A steady state solu- 
tion of Eqs. I|B1|I - P3|I is given by 



Qm,o{t) — 



^ Ki{m) $3(^) + K2{m) ^^iuj) 



(B4) 



im,2 



it) - - 



2/0 

15 |7m 



■<^2{t0) 



R 

-^Kiim) - 13/^2 (m) 
^p 



]m,2n{t) — 



(B5) 
(B6) 



where 71 > 1 in Eqs. ljB6|) . The ratio K^,{m) 
K2 (iTi) / Ki (to) is given by 



i^.(TO)= i^'"(^°^-^ 



riErn{ro) - 1 

(2 + m)(l-rb)flQ 
2mLo 



(B7) 



where 

E,n{ro) = exp 

and — Ri,/Rq. 

The function Cn^^{X) entering in Eq. pSf) has the fol- 
lowing properties: 



(1 - X^) Cf/^iX) C^i'^iX) dX 



{n+l){n + 2) 



(B8) 



n-h3/2 

and cf^iX) = 1, C^^'^iX) = {3/2){5X^ -1). Note that 
due to the condition IjBSp . the function Cq^^{X) only 
contributes to the total angular momentum of the Sun. 
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